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ABSTRACT 

Understanding the thermal equihbrium (stabiUty) curve may offer insights into the nature 
of the warm absorbers often found in active galactic nuclei. Its shape is determined by factors 
like the spectrum of the ionizing continuum and the chemical composition of the gas. We find 
that the stability curves obtained under the same set of the above mentioned physical factors, 
but using recently derived dielectronic recombination rates, give significantly different results, 
especially in the regions corresponding to warm absorbers, leading to different physical pre- 
dictions. Using the current rates we find a larger probability of having thermally stable warm 
absorber at 10^ K than previous predictions and also a greater possibility for its multiphase 
nature, the results obtained with the current dielectronic recombination rate coefficients are 
more reliable because the warm absorber models along the stability curve have computed co- 
efficient values, whereas previous calculations rehed on guessed averages for the same due to 
lack of available data. 
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1 INTRODUCTION 

Warm Absorbers are highly photoionised gas found along our line 
of sight to the continuum source of active galactic nuclei (AGN). 
■ Their signatures are a wealth of absorption lines and edges from 
highly ionized species, notably OVII, OVIII, FeXVII, NeX, CV 
and CVI in the soft X-ray (0.3-1.5 keV) spectra. The typical col- 
umn density observed for the gas is A'^e ~ 10^^*^ cm~^ and the 
temperature T is estimated to be a few times 10^ K. For many ob- 
jects the warm absorber exists as a multiphase absorbing medium 
with all the phases in near pressure equilibrium (Morales, Fabian & 
Reynolds, 2000; CoUinge et al. 2001 ; Kaastra et al. 2002; Krongold 
et al. 2003; Netzer et al. 2003; Ki-ongold et al. 2005; Ashton et al. 
2006). 

Any stable photoionised gas will lie on the thermal equilib- 
rium curve or 'stability' curve, in the temperature - pressure phase 
space, where heating balances cooling; this curve is often used to 
study the multiphase nature of the photoionised gas. The equilib- 
rium depends on the shape of the ionizing continuum and the chem- 
ical abundance of the gas (&olik, McKee & Tarter, 1981; Krolik 
& Kriss, 2001; Reynolds & Fabian, 1995 and Komossa & Mathur, 
2001). We are investigating this dependences in details and will re- 
port elsewhere (Chakravorty et al. , to be submitted). 
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Over the past two decades the estimates of dielectronic re- 
combination rate coefficients have improved. In this letter we show 
that these have affected the stability curves significantly which may 
lead to quite different physical models for the warm absorber gas. 
We conclude with a caution on the reliability of earlier results. 



2 THE STABILITY CURVE 

As is customary, we consider an optically thin, plane parallel slab 
of gas being photoionised by the central source of the AGN. In ion- 
ization equilibrium, photoionisation is balanced by recombinations. 
Thermal equilibrium is achieved when photoionisation and Comp- 
ton heating are balanced by collisional cooling, recombination, 
line-excitation, bremsstrahlung and Compton cooling. The condi- 
tions under which these equilibriums are achieved depend on the 
shape of the continuum, metallicity of the gas, the density, the col- 
umn density, and the ratio of the ionizing photon flux to the gas den- 
sity. Following the convention of Tarter, Tucker & Salpeter (1969) 
we specify this ratio through an ionization parameter ^ — L/ nFI? , 
where L is the luminosity of the source and n the number density 
of gas at a distance R from the center of the AGN. 

We consider a sequence of models for a range of ioniza- 
tion parameters ^, for optically thin gas having constant density 
of lO" cm~^, being irradiated by a power-law ionizing continuum 
with photon index V — 1.8, so that f(i^) ~ and extend- 

ing from 13.6 eV to 40 keV. The chemical composition of the gas 
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Figure 1. Stability curves generated by versions C07 and C84 of CLOUDY, 
for an optically thin gas with constant density 10^ cm~^, being ionized by a 
power-law continuum with photon index V = 1.8, extending from 13.6 eV 
to 40 keV. The gas, in both cases, is assumed to have old Solar abundance 
(see text). Regions of the plane where heating or cooling dominates are in- 
dicated. The curves are seen to differ significantly in the temperature range 
of 4.2 < logT < 7.0. Multiphase regions on the solid and dashed curve 
have been highlighted. 

(referred to as old Solar abundance) is is due to Grevesse and An- 
ders (1989) with extensions by Grevesse and Noels (1993). These 
parameter values (hereafter called the standard set of parameters) 
have been considered because that allows us to compare our results 
easily with earlier work. 

Using version C07.02.01^ (hereafter C07) of the photoionisa- 
tion code CLOUDY (see Ferland et al. 1998 for a description), we 
plot the thermal stability curve which is shown as the solid curve 
(labeled C07) in Figure 1. ^/T is proportional to prad/p, where 
Prad is the radiation pressure and p the gas pressure; so an isobaric 
perturbation of a system in equilibrium is represented by a vertical 
displacement from the curve, and only changes the temperature. If 
the system being perturbed is located on a part of the curve with 
positive slope, which covers most of the curve, then a perturbation 
corresponding to an increase in temperature leads to cooling, while 
a decrease in temperature leads to heating of the gas. Such a gas is 
therefore thermally stable. But if the system is located on one of the 
few parts of the curve with negative slope, then it is thermally un- 
stable because isobaric perturbations will lead to runaway heating 
or cooling. 



3 COMPARISON WITH PREVIOUS WORK 

Reynolds & Fabian (1995, hereafter RF95), have studied the sta- 
bility curve for the warm absorber in MCG-6-30-15, using the 
standard set of parameters given above. We have reproduced their 
curve, as given in Figure 3 of RF95, using version C84. 12a (here- 
after C84) of CLOUDY, which was the stable version between 
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Version ^5 Vpha^es Alog(g/T) AMlog(C/T) 

~ IQS K ~ 10^ K 

C84 45 2 0.05 0.47 0.05 

COT 74 2 0.22 0.46 0.07 



Table 1. Parameters for the warm absorber obtained using versions C84 
and C07 of CLOUDY. The second column gives the value of the ionization 
parameter of the ~ 10® K stable warm absorber phase, the third column the 
number of discrete phases in the multiphase medium, the fourth and fifth 
columns the range of log(^/T) respectively for the ~ 10® K and ~ 10*^ K 
stable warm absorber phases and the sixth column the range of log(5/T) 
over which a multiphase medium is obtained. 



1993 and 1996; this is shown as the dashed curve in Figure 1. 
The C84 and CO? curves match at high temperatures in the range 
T > 10^'^ K, where Compton heating and cooling dominate. They 
also agree at low temperatures T < lO'*'® K. However, the curves 
are significantly different in the intermediate temperature range 
10*'® < T < 10^'^ K, which is the region of interest for the warm 
absorbers and where recombination and line excitation are domi- 
nant cooling mechanisms. 

A detailed comparison between the stable phases for warm 
absorbers predicted by the two different versions of CLOUDY is 
done in Table 1. The second column of the table shows that for an 
absorber at T ~ 10® K, C84 predicts ^5 ~ 45, as compared to 
^5 ~ 74 obtained from the C07 stability curve, where ^5 is de- 
fined to be the ionization parameter corresponding to the middle 
of the 10® K stable warm absorber phase. Both C84 and CO? pre- 
dict two discrete phases of warm absorber at ~ 10® and ~ 10® '' K 
which are in pressure equilibrium with each other and have been 
highlighted in Figure 1. The CO? curve continues to have stable 
thermal states at ~ 10*' K which is not true in the C84 case. The 
range of log(f /T), over which the warm absorber exists are given 
in the fourth and fifth columns of respectively for the low and high 
temperature states. The extent of the 10® K phase in log(^/r) is 
about four times larger in CO?, predicting greater probability of 
finding 10'' K warm absorbers. In the sixth column we have com- 
pared the range Am log(^/r) where the warm absorber exhibits 
multiple phases. CO? predicts a 40% larger range and hence greater 
possibility of a multiphase warm absorber. 

In order to isolate the atomic physics underlying the change 
in the stability curve, we have plotted the fractional variation of 
temperature AT/Tco7 ~ (Tc84 — Tco7)/Tco7, from one version 
to another, against log ^ in the top panel of Figure 2. We see that 
AT > for a major part of the range 1.0 < log^ < 4.5, so 
that the gas is predicted to be cooler by CO?. The cooling fractions 
(AC) of the major cooling agents and the heating fractions (AH) 
contributed by the principal heating agents using CO? are plotted 
against log ^ respectively in the middle and bottom panel of Fig- 
ure 2. The ions that contribute significantly where AT/Tcot 
0.5, are He^^ and high-ionization species of silicon (-1-IO and +11) 
and iron (-1-21, +22 and -1-23). In the same log ^ range, the principal 
heating agents are highly ionized species of oxygen (-1-6 and -I-?) 
and iron (-1-I? to -1-25). 

To identify the ions which are responsible for AT/Tco7 ^ 
50%, we compare their column densities predicted by C84 and 
CO?. In the previous photoionisation calculations in this paper it 
was sufficient to use one zone models for optically thin gas. How- 
ever, to calculate and compare the column densities of the ions over 



DR & warm gas in AGN 3 




log ^ 



Figure 2. Top : The fractional change of temperature between the two versions C84 and C07 of CLOUDY as a function of the ionization parameter. Middle : 
The cooling fraction (AC) for the major cooling agents as a function of log using C07. He+^ and highly ionized species of siUcon and iron play important 
roles. Bottom : The heating fraction (AH) contributed by the principle heating agents as a function of the ionization parameter, obtained using C07. Highly 
ionized species of oxygen like 0+^, 0+^, and various ions of iron hke Fe+^ to Fe+i" and Fe+l'^ to Fe+^s are the key ingredients that are responsible for 
heating the gas in the range 1.0 < log ^ < 4.5. 



the range 1.0 < log^ < 4.5, we chose to specify the gas to have 
total hydrogen column density TVh ~ 10^^ cm~^, which is typical 
for warm absorbers. The column densities are plotted in Figure 3 
with solid lines for CO? and dashed lines for C84. It is seen that 
the column densities of the major coolants changed significantly. 
The cooling agents are among the ions for which dielectronic re- 
combination rate coefficients (hereafter DRRC) have been updated 
to the references below, as will be discussed in detail in Section 4. 
Thus the enhanced cooling in CO? due to the change in DRRC is 
the cause of the shift in the stability curves. 



4 CHANGES IN DIELECTRONIC RECOMBINATION 
RATE COEFFICIENTS 

The evolution of the thermal phases from the nebular temperatures 
of lO'* K to the coronal temperatures of lO" K depends sensitively 
on the detailed atomic physics of the various elements which con- 
tribute to photoelectric heating as well as to cooling due to recom- 
bination and coUisionally excited lines. In Table 2, we have com- 
pared the total recombination rates (dielectronic -I- radiative) for the 
significant cooling agents as predicted by CO? and C84. The val- 
ues for log T and log ^ at which the comparisons have been made 
are given in columns 2 and 3 respectively. For all the ions the total 



Ion 



logT log 5 



Recombination rates (s ^) 



C07 



C84 



He+i 5.34 



2.00 



1.66 



Si+iO 


1.36x10-1 


2.26 xlO-2 




8.85 xlO-2 


2.50 xlO-2 


Fe+21 


9.06x10^2 


2.20 xlO-2 


Fe+22 


8.36x10-2 


2.37 xlO-2 


Fe+23 


5.77 xlO-2 


2.28 xlO-2 



Table 2. Recombination rates of the dominant cooling agents obtained us- 
ing C84 and C07. The ions are given in column 1. The temperature and 
ionization parameter at which the recombination rates have been obtained 
are noted in columns 2 and 3 respectively. The fourth and the fifth columns 
respectively give the value of the recombination rates predicted by C07 and 
C84. 



recombination rates are significantly higher in CO? than in C84 as 
shown by columns 4 and 5. Referring to Figure 2, we see that the 
differences in total recombination rates are relatively larger for ions 
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Figure 3. The column densities for the same ions as in the middle and bottom panel of Figure 2. The solid lines are predicted by C07 and the dashed lines 
are for C84. The difference in the column densities are more pronounced for the ions which are significant cooling agents, with 51+^^^ showing as much as a 
difference of about two orders of magnitude at log g ~ 2.0 where it acts as the dominant cooling ion. 



like Si^^°, Si^^^, Fe"'"'^^ and Fe"'"'^^ which are significant coohng 
agents for log ^ ~ 2.0 and log ^ ~ 3.2, corresponding to which we 
have maximum difference in predicted equilibrium temperatures 
Tco7 and Ten ■ 

In the warm absorber temperature range 10^ < T < 10^ K, 
dielectronic recombination dominates over radiative recombination 
for many ions (Osterbrock & Ferland, 2006). Unhke the radiative 
recombination rate coefficients, the DRRC have undergone signif- 
icant changes over the last decade. C84 used DRRC from Nuss- 
baumer & Storey (1983;1984; 1986; and 1987) and Amaud & Ray- 
mond (1992). In COT we have taken the DRRC for the isoelectronic 
sequences of lithium, beryllium, boron, carbon, nitrogen, oxygen 
and iron-like ions respectively from Colgan et al. (2004), Colgan 
et al. (2003), Altun et al. (2004), Zatsarinny et al. (2004a), Mit- 
nik & Badnell (2004), Zatsarinny et al. (2003), and Gu (2003). The 
DRRC for Ne to Na like ions and Na to Mg like ions are taken from 
Zatsarinny et al. (2004b) and Gu (2004). The CO? DRRC for any 
given ion is usually substantially larger than the C84 DRRC when 
the temperature is much lower than the ionization potential of the 
ion. The significant cooling agents in C07 are among these ions. 
This indicates that the updated DRRC database in C07 is the cause 
of the changes in the stability curve. 

It is seen from Table 2 that the increase in recombination rates 
from C84 to C07 is larger for the lower ionization species. The rea- 
son for the large increase, in general, in the low-temperature DRRC 



(referenced above) is the explicit inclusion of the contribution 
from low-lying (in energy) level-resolved dielectronic recombina- 
tion resonances which have been accurately positioned by reference 
to the observed core energies. Unlike high-temperature dielectronic 
recombination, where the full Rydberg series contributes, low- 
temperature dielectronic recombination is not amenable to simple 
scaling or empirical formulae or guesses. Zatsarinny et al (2003) 
illustrate the radical and erratic changes in the low-temperature 
DRRC on moving between adjacent ions of the same isoelectronic 
sequence and, even, on changing from a term-resolved picture to 
level-resolved one for the same ion. 

The DRRC database is still not complete, specially for the 
lower ionization states. For ions which do not have computed 
DRRC values, C07 uses a solution, as suggested by Ali et al. 
(1991): for any given kinetic temperature, ions that lack data are 
given DRRC values that are the averages of all ions with the same 
charge. The advantage of this method is that the assumed rates are 
within the range of existing published rates at the given kinetic 
temperature and hence cannot be drastically off. However we have 
checked whether the points in the C07 stability curves have com- 
puted DRRC values or guessed average values, concentrating on 
the parts of the curve which have multiphase solutions for the warm 
absorber (10^ < T < 10*^ K) and are different from the C84 sta- 
bility curve. We find that all the ions which act as major cooling 
agents for each of these points in the stability curve have reliable 
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computed DRRC values. Thus, the new data base provides a more 
robust measurement of the various physical parameters involved 
in studying the thermal and ionization equilibrium of photoionised 
gas. 



5 CONCLUSION 

We have shown that stability curves for warm absorbers in AGN 
generated by two versions of the ionization code CLOUDY, C84 
and C07, for the same physical conditions, are substantially differ- 
ent in shape, leading to different conclusions regarding the nature 
of the warm absorber. The differences in the results of the photoion- 
isation calculations arise due to major changes in the dielectronic 
recombination rate coefficient (DRRC) data bases which have taken 
place over the last decade. The modern version C07 includes reli- 
able computed DRRC values for many more ions than C84 for the 
entire part of the stability curve relevant for warm absorbers and 
does not rely on guessed average values. Thus, the physical nature 
of the warm absorber predicted by modem calculations are more 
reliable. We suggest that past calculations for photoionised plasma 
in thermal equilibrium should be reconsidered, taking into account 
these changes. Since, any other atomic physics code (eg. XSTAR), 
in popular use for warm absorber studies, is also likely to be af- 
fected by the changes in the DRRC data base, caution should be 
exercised while using or comparing results from them as well. 
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